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Abstract. Existence of different kinds of synchronizations, namely anticipatory, complete 
and lag type synchronizations (both exact and approximate), are shown to be possible in time- 
delay coupled piecewise linear systems. We deduce stability condition for synchronization of 
such unidirectionally coupled systems following Krasovskii-Lyapunov theory. Transition from 
anticipatory to lag synchronization via complete synchronization as a function of coupling delay 
is discussed. The existence of exact synchronization is preceded by a region of approximate 
synchronization from desynchronized state as a function of a system parameter, whose value 
determines the stability condition for synchronization. The results are corroborated by the 
nature of similarity functions. A new type of oscillating synchronization that oscillates between 
anticipatory, complete and lag synchronization, is identified as a consequence of delay time 
modulation with suitable stability condition. 



1. Introduction 

Synchronization of coupled chaotic systems is a field of growing interest in view of its potential 
applications in diverse area of research activities such as secure communication, cryptography, 
controlling, long term prediction of chaotic systems, etc. 0121 01^0 El- Since the identification 
of chaotic synchronization different kinds of them have been proposed in interacting chaotic 
systems both theoretically and experimentally: complete synchronization [7| [S], generalized 
synchronization j^l phase synchronization I12j . lag synchronization jl3l [H] and 

anticipatory svnchronizationl lSI I16j . 

One of the most important applications of chaos synchronization is secure communication. It 
has been shown that in secure communication based on simple low dimensional chaotic systems 
with only one positive Lyapunov exponent, the hidden message can be unmasked by dynamical 
reconstruction of the chaotic signal using nonlinear dynamical forecasting methods or by using 
some simple return maps |171 I18j . One way to overcome this problem is to consider chaos 
synchronization in higher dimensional systems having multiple positive Lyapunov exponents, 
based on the consideration that increased randomness and unpredictability of the hyperhcaotic 
signals will make it more difficult to extract the masked message. Recently, chaotic time-delay 
systems have been suggested as good candidates for secure communication '19, as the time- 
delay systems exhibit the intriguing characteristics of increase in the embedding dimension and 
the number of positive Lyapunov exponents with the time-delay inspite of its small number 



of physical variables. Therefore the study of chaos synchronization in time-delay systems is of 
considerable practical significance. 

Recently, we have shown that even a single scalar delay equation with piecewise linear 
function can exhibit hyperchaotic behavior even for small values of time-delay and hence it 
is of interest to consider chaos synchronization in such scalar piecewise linear delay differential 
equations with appropriate coupling between them. For the present study, we will consider 
chaos synchronization between two such time-delay systems with unidirectional coupling and 
having two different time-delays: one in the individual systems, namely, feedback delay ti and 
the other in the coupling term T2 ■ We have arrived at the stability condition for synchronization 
following Krasovskii-Lyapunov theory, which shows that the stability condition is independent 
of the delay times, and demonstrate that there exist transitions between three different kinds 
of synchronization, namely, anticipatory, complete and lag synchronizations by simply tuning 
the second time-delay parameter in the coupling, for a fixed set of other system parameters 
satisfying appropriate stability condition. The results have been corroborated by the nature 
of similarity functions. To enhance the security, delay time modulation is introduced at the 
feedback delay and the stability of the synchronization manifold with delay time modulation is 
calculated. We have also identified a new type of oscillating synchronization as a consequence 
of delay time modulation. 



2. Scalar piecewise linear delay system 

We consider the following first order delay differential equation introduced by Lu and He |2j 
and discussed in detail in references [201 [221 122] , 

x{t) = -ax{t) + hf{x{t-T)), (1) 

where a and b are parameters, r is the constant time-delay and / is an odd piecewise linear 
function defined by the relation 
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This system exhibits hyperchaotic behavior for the parameter values a = 1.0,6 = 1.2 and r = 5 
and the hyperchaotic nature was confirmed by the existence of multiple positive Lyapunov 
exponents. The first ten maximal Lyapunov exponents as a function of r for the chosen 
parameter values is shown in Fig. ^ In our present study, we have included three linear parts 
for the function f{x) in Eq. [21 One may also carry out the studies with two linear parts alone. 
However, the dynamics in the three linear part case is richer and so we confine our study to this 
case alone here. 



3. Stability condition for synchronization 

Now let us consider the following unidirectionally coupled drive xi{t) and response X2{t) systems 
with two different time-delays ri and T2 as feedback and coupling time-delays, respectively, 

XI (t) = -axi(t) + 6i/(xi(t - n)), (3) 



X2(t) = -ax2{t) + b2f{x2{t - n)) + bsfixiit - T2)), 
where 61, 62 and 63 are constants, a > 0, and /(x) is of the same form as in Eq. Q. 



(4) 




Now we can deduce the stability condition for synchronization of the two time-delay systems 
Eqs. © and © in the presence of the delay coupling h^f{xi{t — T2)). The time evolution of the 
difference system with the state variable A = xir — X2 , where xir = xi{t — t),t = T2 — ti, can 
be written for small values of A by using the evolution Eqs. (3) as 

A = -aA + (62 + 63 - h)fixiit - T2)) + b2f'ixiit - T2))Ar,, (5) 

In order to study the stability of the synchronization manifold, we choose the parametric 
condition, 

bi = b2 + h, (6) 
so that the evolution equation for the difference system A becomes 

A = -aA + 62/'(xi(t-r2))A,,. (7) 

The synchronization manifold is locally attracting if the origin of this equation is stable. 
Following Krasovskii-Lyapunov functional approach [24| I25j . we define a positive definite 
Lyapunov functional of the form 

vit) = ^a2 + /i r A2(t + e)de, (8) 

where n is an arbitrary positive parameter, /i > 0. Note that V{t) approaches zero as A ^ 0. 

To estimate a sufficient condition for the stability of the solution A = 0, we require the 
derivative of the functional V{t) along the trajectory of Eq. 0, 

^ = -aA2 + b2f'{xi {t - T2))AAr, + fiA^ - fiA^ , (9) 

to be negative. After simple algebra, one can find the sufficient condition for asymptotic stability 
as 

a>\b2fixiit-T2))\ (10) 
along with the condition ® on the parameters 61, 62 and 63. 



Now from the form of the piecewise hnear function /(x) given by Eq. (2), we have, 




1.5, 0.8<|xi|<| 
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(11) 



Consequently the stabihty condition (|lflj) becomes a > 1.5|52| > I ^2 1 along with the parametric 
restriction bi =62 + ^3- 

Thus one can take a > I62I as a less stringent condition for (fTfl|) to be valid, while 



as the most general condition specified by ()1U() for asymptotic stability of the synchronized state 
A = 0. The condition (|12() indeed corresponds to the stability condition for exact anticipatory, 
identical as well as lag synchronizations for suitable values of the coupling delay T2- It may also 
be noted that the stability condition (fT^ is independent of the both the delay parameters ti 
and T2. 

4. Anticipatory synchronization 

To start with, we demonstrate that there exists a region of anticipatory synchronization when 
T2 < Ti at the synchronization manifold. For this purpose, we have fixed the value of feedback 
delay ri at ri = 0.25 while the other parameters are fixed as a = —0.16, 61 = 0.2, ^2 = 0.1, 63 = 
0.1 and the coupling delay T2 is considered as a control parameter. For all values of T2 < ti, 
we find from the evolution equation the difference variable A = xi{t — t) — X2{t) with 
T = T2 — T1 < approaches zero asymptotically, so that xi(t + r) = X2(i), thereby demonstrating 
the existence of anticipatory synchronization with anticipating time as r. The existence of 
anticipatory synchronization is characterized by the nature of similarity functions for different 
values of 62, the parameter whose value determines the stability of the synchronization manifold 
as may be seen from Eqn. (|12|) . 

Now we use the notion of similarity function to characterize anticipatory synchronization 
which was originally introduced by Rosenblum et al. to characterize lag synchronization. 
Similarity function is defined as a time-averaged difference between the variables xi and X2 (with 
mean values being subtracted) taken with the positive time shift r for lag synchronization and 
negative time shift — r for anticipatory synchronization. 



where {x) means time average over the variable x. If the minimum value of Sa{T) reaches 
zero, that is Sa^r) = 0, then there exists a time shift — r between the two signals xi{t) and 
X2{t) such that X2{t) = xi(t + r), demonstrating the existence of anticipatory synchronization 
between the drive xi and the response X2 signals. Fig. |21 shows the similarity function 5a (r) 
as a function of the coupling delay T2 for four different values of 621 the parameter whose value 
determines the stability condition given by Eq. H12|l . while satisfying the parametric condition 
bi = b2 + 63. Curves 1 and 2 are plotted for the values of 62 = 0.18(> a = 0.16 > a/1.5) 
and 62 = 0.16(= a > a/1.5), respectively, where the minimum values of Sa(r) is found to be 
greater than zero, indicating that there is no exact time shift between the two signals xi{t) 
and X2{t). Note that in the both cases the stringent stability condition p2j) and the less 
stringent condition a > I62I are violated. Curve 3 corresponds to the value of ^2 = 0.15 ( 
< a > a/1.5), where the minimum value of 5'a(r) is almost zero, but not exactly zero (as may 
be seen in the inset of Fig. [21 inwhich case only the less stringent condition a > I62I is satisfied 
while the stringent stability condition is ruled out), indicating an approximate anticipatory 
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Figure 2. Similarity function 
5a (r) for different values of b2, the 
other system parameters are a = 
0.16, 61 = 0.2 and n = 25.0. 
(Curve 1: 62 = 0.18,63 = 0.02, 
Curve 2: 62 = 0.16,63 = 0.04, 
20 22 24 26 28 30 Curve 3: 62 = 0.15,63 = 0.05 and 

Curve 4: 62 = 0.1,63 = 0.1). 

synchronization xi{t) ~ X2{t — t). On the other hand the curve 4 is plotted for the value 
of 62 = 0.1(< a/1.5), satisfying the general stability criterion, Eq. H12|) . It shows that the 
minimum of Sair) = 0, thereby indicating that there exists an exact time shift between the two 
signals demonstrating anticipatory synchronization. The anticipating time is found to be equal 
to the difference between the coupling and feedback delay times, that is, r = T2 — ti. Note 
that (t) = for all values of T2 < ri, indicating anticipatory synchronization for a range of 
delay coupling. A further significance is that the anticipating time t = T2 — ti is an adjustable 
quantity as long as T2 < ri, which can be tuned suitably to satisfy experimental situations. 
Thus anticipating of any future chaotic state can be made possible. We also note that the above 
choice of parameters (Curve 4) in Fig. |21is only a specific example exhibiting exact anticipatory 
synchronization. Any other choice of parameters satisfying the condition 61 = 62 + 63 subject to 
a > 1.5 1 62 1 is good enough to obtain the above phenomenon. 

The existence of exact synchronization is preceded by a region of approximate synchronization 
from desynchronized state as a function of the parameter 62, whose value determines 
the stability condition for synchronization. We have also found !231that the existence of 
approximate synchronization from desynchronized state is characterized by transition from on- 
off intermittency to periodic structure in probability distribution of laminar phase as shown by 
Zhan et ai[Tl]. 

5. Complete synchronization 

Now, we show that when the value of coupling delay T2 equals the value of feedback delay ri, 
that is, T2 = Ti, there exists complete synchronization. We find from the evolution equation ((7|, 
the difference variable A = xi{t — T) — X2{t) with t = T2 — T1 = approaches zero asymptotically, 
so that xi{t) = X2{t), thereby demonstrating the existence of complete synchronization between 
the drive xi{t) and the response X2{t). Here also, we have identified that the emergence of 
approximate complete synchronization is associated with the transition form on-off intermittency 
to periodic structure in the laminar phase distribution as a function of the parameter 62 . 

6. Lag synchronization 

Next, we show that for the value of coupling delay T2 greater than feedback delay ri, the system 
exhibits lag synchronization. For T2 > ti, we find from the evolution equation the difference 
variable A = xi{t — t) — X2{t) with t = T2 — ti > approaches zero asymptotically, so that 
xi(t — r) = X2{t), thereby demonstrating the existence of lag synchronization with lag time equal 





to T = T2 — Ti. The existence of lag synchronization is characterized by the similarity function 
Si{t) defined earlier in Eqn. (|TH|) . Fig. |31 shows the similarity function Si{t) Vs coupling delay 
T2 for four different values of 62. Curves 1 and 2 show the similarity function Si{t) for the 
values of 62 = 0.18 and 0.16, respectively. The minimum of similarity function Si{t) occurs 
for values of Si{t) > 0, where the minimum does not occur at zero and hence there is a lack 
of exact lag time between the drive and response signals indicating asynchronization. Curve 3 
corresponds to the value of 62 = 0.15 (which is less than a but greater than a/1.5), where the 
minimum values of Si{t) is almost zero, but not exactly zero (as may be seen in the inset of 
Fig. 121), so that xi{t) ^ X2{t + t). However for the value of 63 = 0.1, for which the general 
condition (|12j) is satisfied, the minimum of similarity function becomes exactly zero (Curve 4) 
indicating that there is an exact time shift (Fig. 13]) between drive and response signals xi{t) and 
X2{t), respectively, confirming the occurrence of lag synchronization. Again the above choice of 
parameters (Curve 4) in Fig. |21 is only a specific example exhibiting exact lag synchronization. 
Any other choice of parameters satisfying the condition © subject to a > 1.5|62| is good enough 
to obtain the above phenomenon. 

We have also confirmed that as in the case of anticipatory synchronization, when the 
parameter 62 varies, the onset of exact lag synchronization is preceded by a region of approximate 
lag synchronization, which is characterized by a transition from on-off intermittency of the 
desynchronized state to a periodic structure in the laminar phase distribution. 

7. Stability condition for synchronization with delay time modulation 

To improve the security of the above type of synchronized systems, we have extended our studies 
by choosing time-delay as a function of time. Recently, the concept of delay time modulation 
was introduced by Kye et al [5^1 HZ] and they have shown that reconstruction of phase space of 
the time-delay systems is hardly possible as a consequence of delay time modulation. 

Now, we will introduce the delay time modulation in the systems Q and (^J at the feedback 
delays ri (t) in the form 

n(i) = To + TaSm{uit), (14) 

where tq is the zero frequency component, is the amplitude and lv/2tt is the frequency of the 
modulation. Here the coupling delay T2 is kept constant. 

The stability condition for synchronization of the systems (jHl) and (HJ with delay time 
modulation can be obtained as (by following the same procedure as for the case of constant 
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Figure 4. Oscillating synchronization from lag to anticipatory synchronization via complete 

synchronization for tq = 100, = 90 and to = 10~^. Drive xi{t) is represented by and 

response X2{t) by — • — . 



feedback delay given in section 3), 
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along with the condition © on the parameters 61,62 and 63. Once again as in the previous 



case, from the form of piecewise function f{x), one can take a > 
condition for (fT3)) to be valid, while 
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as the most general condition specified by (fT^ for asymptotic stability of the synchronized state. 



8. Oscillatory synchronization as a consequence of delay time modulation 

We have fixed the value of tq = 100, = 90 and lu = 10~^ with the same value of a 
and b as previously studied. For the value of 62 = 0.01 the general stability condition (|15j) 
is satisfied. For the chosen values of tq and r^, one can find that ti oscillates between 
(ri(t) = To + Ta sin(a;t) = 100 it 90) 10 and 190. Hence for any value of coupling delay T2 less 
than 10, T2 is always less than Ti{t), one can find from the evolution equation the difference 
variable A = xi{t — T{t)) — X2{t) with T(t) = T2 — Ti(t) approaches zero asymptotically, so 
that xi{t + T{t)) = X2{t), implying the existence of anticipatory synchronization with the time 
dependent anticipating time T(t) = T2 — ri(t). Similarly for any value of coupling delay T2 greater 
than 190, that is for T2 always greater than Ti(t), we have xi{t — T(t)) = X2{t), demonstrating 
the existence of lag synchronization. On the other hand, when the coupling delay T2 is chosen in 
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Figure 5. Oscillating synchronization from anticipatory to lag synchronization via complete 

synchronization for tq = 100, Ta = 90 and io = 10~^. Drive xi{t) is represented by and 

response X2{t) by — • — . 

between 10 and 190, one is able to find that the feedback delay ti (t) is less than T2 for some time 
(inwhich case r(t) = T2 — Ti(t) > 0, so that there exists lag synchronization xi{t — r(t)) = X2{t)) 
and as Ti{t) increases eventually, it approaches X2{t) at a certain time, where r(t) = T2—Ti{t) = 0, 
so that xi{t) = X2{t), that is there exists complete synchronization for a specific value of 
time. As Ti(t) exceeds the value of r2, now r(t) becomes negative, T(t) = T2 — Ti{t) < 
with xi{t + r(t)) = X2(t)), there exists anticipatory synchronization. Therefore as time evolves 
there is oscillation between lag, complete and anticipatory synchronization with time dependent 
anticipating and lag times. Fig. 0] shows the evolution of the drive xi(t) and the response 
X2{t) at the transition between lag to anticipatory synchronization via complete synchronization, 
whereas Fig.lSJshows the evolution of the drive xi{t) and the response X2{t) at the next transition 
between anticipatory to lag via complete synchronization. Thus as a consequence of delay time 
modulation there exists a new type of oscillating synchronization that oscillates between lag, 
complete and anticipatory synchronization with varying anticipating and lag times. 

9. Summary and conclusion 

In summary, we have obtained the stability condition for synchronization in a system of two 
coupled piecewise linear differential equation following Krasovskii-Lyapunov functional approach 
for both the case of constant and time dependent delay. We have demonstrated that there exists 
regimes of anticipatory, complete and lag synchronizations with coupling delay T2 as control 
parameter for constant delay case with anticipating and lag time as t = T2 — ti. The exact 
synchronization regimes are preceded by a regime of approximate synchronization from the 
desynchronized state as a function of parameter b2, which determines the stability condition. 
The existence of anticipatory and lag synchronizations are characterized by the similarity 
functions. In the case of time dependent delay, we have demonstrated that there exists regimes 
of anticipatory and lag synchronization with time dependent anticipating and lag times for 
suitable values of coupling delay T2. We have also shown the existence of new type of oscillating 
synchronization that oscillates between lag, complete and anticipatory synchronizations as a 



consequence of delay time modulation, when the value of coupling delay T2 is chosen in between 
Ti = To — Ta and Ti = To + Ta with varying anticipating and lag times. We believe that this new 
type of oscillating synchronization with time dependent anticipating and lag times will enhance 
the security of secure communication and now we are implementing this experimentally through 
electronic circuits. 
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